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Abstract 

This is the second of two papers where we study the asymptotics of the gen- 
erahzed Nariai solutions and its relation to the cosmic no-hair conjecture. In the 
first paper, the author suggested that according to the cosmic no-hair conjecture, 
the Nariai solutions are non-generic among general solutions of Einstein's field equa- 
tions in vacuum with a positive cosmological constant. We checked that this is true 
within the class of spatially homogeneous solutions. In this paper now, we continue 
these investigations within the spatially inhomogeneous Gowdy case. On the one 
hand, we are motivated to understand the fundamental question of cosmic no-hair 
and its dynamical realization in more general classes than the spatially homogeneous 
case. On the other hand, the results of the first paper suggest that the instability of 
the Nariai solutions can be exploited to construct and analyze physically interesting 
cosmological black hole solutions in the Gowdy class, consistent with certain claims 
by Bousso in the spherically symmetric case. However, in contrast to that, we find 
that it is not possible to construct cosmological black hole solutions by means of 
small Gowdy symmetric perturbations of the Nariai solutions and that the dynamics 
shows a certain new critical behavior. For our investigations, we use the numerical 
techniques based on spectral methods which we introduced in a previous publication. 

1 Introduction 

In this work, we are interested in a particular consequence of the cosmic no-hair conjec- 
ture [T6l|T9]. As discussed in our first paper [5], this conjecture suggests that the so-called 
Nariai solutions are non-generic in the class of cosmological solutions of Einstein's field 
equations in vacuum 

Gf,^ + Agf,^ = 0, (1.1) 

with a positive cosmological constant A, due to its extraordinary asymptotics for large 
times. Throughout the paper, a cosmological solution means a globally hyperbolic solu- 
tion of Eq. (jl.ip with compact Cauchy surfaces. 
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In the first paper [5], we analyzed the asymptotics of what we cahed generahzed 
Nariai solutions. A particular solution in this family is the (standard) Nariai solution 
|23[ . Although these solutions are all isometric locally to the standard Nariai solu- 
tion, we were motivated to introduce this family for the following reasons. First, their 
spatially homogeneous perturbations have interesting properties [5j. Second, the ap- 
proach presented in this paper here makes particular use of the non-standard Nariai 
solutions in this family. In the following, we often speak of "Nariai solutions", when 
we mean any generalized Nariai solution. A particular contribution of the first paper 
was a proof of the outstanding fact that the Nariai solutions do not possess smooth 
conformal boundaries, a result closely related to the cosmic no-hair picture as explained 
there. Moreover, we investigated how the expected non-genericity of the Nariai solutions 
is realized dynamically in the spatially homogeneous class of perturbations. In general, 
when we speak of a perturbation of a Nariai solutions, we mean a cosmological solution 
of the fully non-linear Einstein's field equations Eq. (jl.ip whose data, on some Cauchy 
surface, is close to the data on a Cauchy surface of a generalized Nariai solution. By 
"close" we mean that two data sets should deviate not too much with respect to some 
reasonable norm in the initial data space. We show in [5] that an arbitrary small spa- 
tially homogeneous perturbation of any Nariai solution does either not expand at all in, 
say, the future, or it expands in a manner consistent with the cosmic no-hair picture by 
forming a smooth future conformal boundary. 

Certainly, the case of spatially homogeneous perturbations is special and it would 
be interesting to study the instability of the Nariai solutions within more general classes 
of perturbations. Beyond the problem of cosmic no-hair, however, it is a tempting 
possibility to exploit our knowledge about the instability in the homogeneous case in 
order to construct new non-trivial inhomogeneous cosmological black hole solutions. 
Recall from the results in the first paper ^ that in the spatially homogeneous case, the 
sign of the initial expansion H^^^ of the spatial §^-factor of a perturbation of a Nariai 
solution controls whether the spatial S^-factor collapses or expands to the future. Hence, 
one can expect that by making H^^^ spatially dependent on the initial hypersurface, we 
become able to control the spatially local behavior of the perturbations. In particular, 
we should obtain solutions with arbitrary many black hole interiors on the one hand 
and expanding cosmological regions on the other hand. In principle, we are interested 
in studying generic inhomogeneous perturbations of the Nariai solutions without any 
symmetries. However, this is not feasible in practice. A systematic approach would 
be to reduce the symmetry assumptions step by step. The first systematic step is the 
spherically symmetric case. Indeed Bousso in [7J claims that arbitrarily complicated 
spherically symmetric cosmological black hole solutions can be constructed with this 
approach. 

In this paper, we do not consider the spherically symmetric case, but rather proceed 
with the Gowdy symmetric case. Gowdy symmetric solutions with spatial T^-topology 
have been very prominent in the field of mathematical cosmology, and recently, important 
outstanding problems have been tackled rigorously; see the important work in \26\ [27j 
based on a long list of previous references. However, the Gowdy case with spatial 
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S"*^ X S^-topology (and §^), which is the relevant case here, has turned out to be more 
difficult analytically [211 US [2H] • This is one particular motivation for us to proceed 
with this class by means of numerical techniques. In fact, we have developed numerical 
techniques in [6] applicable to the x §^-Gowdy class, which will be used in this work 
here. An alternative numerical approach for Gowdy solutions with §^ x §-^-topology can 
be found in [15j. However, this approach is not applicable directly for the conformal 
field equations based on orthonormal frames which we have decided to work with. In 
any case, our results can be hoped to complement the claims in [7] in a physically and 
technically interesting setting. 

In all what follows, we use the same fundamental conventions and assumptions as in 
the first paper [5]. 

Our paper is organized as follows. In Section [21 we prepare our investigations. After 
a short introduction to Gowdy symmetry on S X §2 in Section \2.1\ we construct those 
Gowdy invariant initial data sets in Section 12.21 which will be used to study the pertur- 
bations of the Nariai solutions later. For the basic properties of these solutions, we refer 
to Section [3] of [5j. Then, Section [2T3l of the paper here is devoted to the discussion of 
certain mean curvature quantities which are analogous to those quantities which control 
the instability of the Nariai solution in the spatially homogeneous case. In Section 12.41 
we present the formulation of Einstein's field equations that we will use for the numerical 
computations, namely the conformal field equations. We briefly describe the unknown 
variables in this formulation and fix our choice of gauge. We conclude with a comment 
about our particular reduction to 1 + 1. We proceed by giving a quick summary of the 
numerical infrastructure in Section 12.51 In Section 12.61 we comment on how to compute 
initial data for the conformal field equations from the data constructed in Section 12.21 
Numerically, this is not completely trivial, because the initial data before are based on 
coordinate components of the metric and hence they run the risk of coordinate singu- 
larities when they are transformed to orthonormal frame based data for the conformal 
field equations. In Section 13.11 we fix and motivate particular initial data sets for the 
later numerical runs. The central part of the paper is Section 13.21 where we show the 
numerical evolutions and interpret the results. In Section [3.31 we present further details 
of a practical nature about the numerical runs and then proceed with an analysis of nu- 
merical errors. The paper is concluded with a summary, a discussion of open problems 
and an outlook in Section [H 

2 Preparations 

2.1 Gowdy symmetry on §^ x 

We quickly introduce important and relevant facts about Gowdy symmetry on x S^. 
General aspects of U(l) x U(l)-symmetric solutions of Einstein's field equations were 
discussed in p[7] for the first time and later reconsidered in \9\. 
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Smooth U(l) X U(l)-invariant metrics on S-*^ x §^ Let us introduce coordinates 
{p,9,(p) on X where p G (0,27r) is the standard parameter on and {9,(j)) are 
standard polar coordinates on S^. The coordinate vector fields dp and d^p generate a 
smooth effective action of the group U(l) x U(l) on x S^. Let us consider a smooth 
Riemannian metric h on §^ xS^ which is invariant under this action. One can parametrize 
h as 

h = e^^dO'^ + R{e^d(t)'^ + 2e^Q# dp + (e^Q^ + e-^)dp^). (2.1) 

In particular, the field dg can be assumed to be orthogonal to the group orbits every- 
where. All functions involved here only depend on 9. For the smoothness of /i, it is 
sufficient and necessary [S] that there are function i?, P and A so that 

R = Rsm9, P = P + lnsm9, X = {P + lnR)/2 + X, (2.2a) 

where R, P, X and Q are smooth function of cos^, and 



Gowdy symmetric spacetimes with spatial §^ x §^-topology Now let us con- 
sider a globally hyperbolic spacetime {M,g) foliated with U(l) x U(l)-invariant Cauchy 
surfaces of topology x S^, in the sense that the first and second fundamental form of 
each surface are invariant under the U(l) x U(l)-action before. It has been shown before 
[21 [To] that this action is orthogonally transitive, i.e. the twist constants 

vanish for spatial topology x (or S^, but not T^) if this spacetime is a solution 
of Eq. (jl.ip . Here all index manipulations are done with the metric g, and e is the 
volume form of g assuming that the spacetime is orientable. See also |29] for more 
details. In general, if these constants vanish, a U(l) x U(l)-invariant spacetime is called 
Gowdy spacetime [17]. Thus, all U(l) x U(l)-invariant solutions of Eq. (jl.ip with spatial 
topology S"*^ X (or S^) are Gowdy solutions. 

It is clear that that there exist coordinate gauges which are inconsistent with the 
assumption, that for all t = corest-hypersurfaces, the first and second fundamental form 
are U(l) x U(l)-invariant and that the Killing vector fields can be identified with the 
coordinate vector fields dp and 5^ for all values of the time coordinate t. Most prominent 
examples of - in this sense - consistent gauges are the so-called areal gauge where one 
sets R = t, and the "conformal time gauge"LJ where the time coordinate t is chosen so 
that 

g = -e^^dt^ + h (2.3) 

^The name "conformal time gauge" must not be confused with the conformal approach described in 
Section [2.41 This name was chosen because the 2-surfaces orthogonal to the group orbits are explicitly 
conformally flat. 
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with h given by Eq. (j2.1|) together with Eqs. (j2.2|) : for more information see for instance 
[lOj . A further "consistent" choice of gauge is the Gauss gauge. Here the time coordinate 
t is chosen so that metric g takes the form 

g = -dt^ + h (2.4) 

with the same h as before. 

Spatial homogeneity and the Nariai case If {M,g) is spatially homogeneous, it 
is in particular Gowdy symmetric, and hence the metric takes the following form in the 
conformal time gauge 

g = Re^{-dt^ + de'^ + sin^ Od^'^) + Re-^dp^. (2.5) 

All functions in this metric only depend on t. In this gauge, the generalized Nariai 
metrics are determined by 

R{t) = ^{t)/K, P{t) = -ln^{t). (2.6) 

2.2 A family of Gowdy symmetric initial data close to Nariai data 

Our aim is now to find Gowdy invariant initial data close to the Nariai solution. These 
will be interpreted as perturbed Nariai data, and the corresponding solutions of the 
field equations as perturbed Nariai solutions. Our data sets must be solutions of the 
constraint equations implied by the vacuum Einstein's field equations with A > on a 
Cauchy surface of x §^-topology. 

Since Gowdy symmetry is generated by the coordinate vector fields dp and c^^, our 
initial value problem reduces to a problem on the domain of the coordinates (t, 8) in 
principle. Let t be the time coordinate of the "conformal time gauge" defined in Eq. ()2.3p . 
The constraints are 

• Hamiltonian constraint: 

1 „,2 1 9P^;2 R'^ 1 ^9 1 9P^9 9\ . R' X' R" RX 

• Momentum constraint 

U 2^^+ 2^ 2ii2 R R ^ R- 

A dot represents a t-derivative and a prime a ^-derivative. We assume that on the 
initial hypersurface, the action of the Gowdy group is of the standard form, which 
implies that the quantities R, P, Q, X must be expressible via Eqs. (|2.2p . When the new 
quantities R, P and A are substituted into the constraint equations, formally singular 
terms arise at the coordinate singularities 9 = 0, vr. In this paper, we do not address 
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the problem of these terms; another future publication will be devoted to such and 
related issues. In the case A = 0, the corresponding problem arises [15j. For A > 0, 
however, Eq. (14) in [15] must be substituted by a more complicated condition due to 
the additional term e^^A. In order to circumvent these problems for the time being as 
in [15], we are satisfied with a particular family of explicit solutions of the constraints 
with the properties above for this paper. However, it is clear that such a family of initial 
data cannot be considered "generic" and hence no strict results about the cosmic no-hair 
conjecture can be expected. Nevertheless, we see the investigations in this paper as a 
promising first step. 

In all of what follows, we assume A = 3 with loss of generality, since it yields the 
simplest expressions and makes all quantities dimensionless. We have mentioned before 
that the functions R, P, Q and A must be smooth functions of cos 9, and A has to become 
zero at 9 = 0, vr. Now, we make a polynomial ansatz for these functions in z = cos 9, 
and solve the constraints for the polynomial coefficients matching these conditions. In 
this way, which requires cumbersome algebra done with Mathematica, we derive the 
following almost explicit family of Gowdy symmetric solutions of the constraints 

R = R^, R! = — , (2.7a) 

P = P^- ^N^^) W 9, P' = ^(SL°^ - S« cos^ 9), (2.7b) 

Q = _^iv« r\~n^)dz, Q' = ^sL'^e-^, (2.7c) 



K J _i K 



A = ^ivW W 9, A' = ^(V3sf - S«) sin2 9. (2.7d) 

The constants P=k , , are determined transcendentally by the initial data param- 
eters (i?^=, K, ^x^) in the following manner. Let us make the abbreviation 

C := ^J^' + (3(S«)2 - 6SL°)S« + 2V3SL°^) + 3(S« - sL°))2. 



Then 



P* =ln 



1 /2^ _ 2C. _ ^^(1)^2 _ ^5.(0))2 + + 1 



4P,k2 \^ 3 ^ 



K 



K V3 

^ ^ Si^i^r - 6S(°)S« + 3Ck(s(°) - S«) + (VSk' + 3S(°) 
5^(2) 



3^2 

The reality conditions on both the root in the definition of C and the logarithm in 
the definition of P* imply restrictions for the choice of the (otherwise free) parameters 
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(i?*, K, S_ , ) which, however, we do not make exphcit now. For the apphcations 
later, we always check that these are satisfied without further notice. Let us suppose 
that all quantities in Eqs. (|2.7|) are well defined. The only non-explicit expression is the 
integral for Q when t!"^ ^ 0. We compute this integral numerically by approximating 
the exponential by its truncated Taylor series. This series converges very quickly and in 
practice, the series can be truncated after a few terms. 

Let us also remark that these data are not polarized in general, i.e. the Killing fields 
cannot be chosen globally orthogonal. 

Now, let us identify spatially homogeneous, and in particular Nariai data in our 
family Eq. ()2.7p . Writing Eq. ()2.5p for the conformal time gauge according to Eq. ()2.3p . 
we see that spatial homogeneity implies S^J^^ = N^"^ = = 0. Using the expressions 
above, the data are hence spatially homogeneous if and only if 

In particular, t!^^ plays the role of an "inhomogeneity parameter". However, our family 
of data does not comprise all spatially homogeneous data, in particular not all Nariai 
data. For any generalized Nariai data, Eq. ()2.6p implies 

= $*/A, K = = -l/\/3, t!-^^ = 0. (2.9) 

In this case, Eq. (j2.8p yields < 1, and thus only generalized Nariai solutions with 
(To < are present in our data. 

2.3 Mean curvatures and the expected instability in the Gowdy class 

In principle, there is no "canonical" quantity which could play the same role for the 
instability of the Nariai solutions in the spatially inhomogeneous case as the quantity 
hI^^ defined in Section 14.11 in the first paper [3] in the spatially homogeneous case. 
First, there is no "canonical" foliation of spacetime, and, second, no "geometrically 
preferred spatial S^-factor". Here, we choose a Gaussian foliation with time coordinate 
t; cf. Eq. ()2.4p . At this stage, we can only hope that at least for small perturbations of 
the Nariai solutions, the "expansion of the coordinate S-^-factor", which we define now, 
plays a similar role here as ; at least these quantities agree for spatially homogeneous 
perturbations. 

On the initial hypersurface, we choose coordinates (p, 0, (j)) as before. By means of 
the Gauss gauge condition, these spatial coordinates are transported to all t = const- 
hypersurfaces S^. On any S^, a "coordinate S^-factor" is then a 2-surface diffeomorphic 
to §^ determined by p = const. Since the metric is invariant under translation along p, all 
such 2-surfaces are isometric at a given t. Similarly, we define "coordinate S^-factors"; 
note that these 1-surfaces are not isometric on a given S^. The expansions (mean 
curvatures) associated with these surfaces are defined as follows for a Gauss gauge. Let 
H be the mean curvature of any t = consi-hypersurface. Let H2 be the projection of the 

7 



(2.8) 



mean curvature vector of a coordinate S'^-factor to dt, this is the quantity we refer to as 
the "expansion of the coordinate §^-factor" . Similarly, we define Hi as the "expansion 
of the coordinate S^-factor". With the expression Eq. ()2.ip for the spatial metric, the 
following formulas hold 

m' + r(p' + 2A'l r' + r(p' + X' 



H = H2 = (2.10a) 

6R 2R 

e^^Qsm^eiqP' + Q'\ 
3H -2H2- Hi = ^, (2.10b) 

where a prime denotes a derivative with respect to Gaussian tim^. 

We expect that, as for the quantity H^^^ in the spatial homogeneous case, the sign 
of the initial value of H2 controls whether the solution collapses or expands locally in 
space and hence plays a particularly important role for the description of the expected 
instability of the Nariai solutions. Thus we write down the expression for the family of 
initial data in Eqs. (j2.7p 



-"2 1 initial ^^3 



3(S« - - 3kC(S« - SL°)) + ^/3k2(S« + cos^ ^ ^ ^^'^^^ 



2.4 Formulation of Einstein's field equations 

Our numerical approach, based on orthonormal frames for the conformal field equations 
is discussed in [6l |3] and is briefly summarized in the following paragraphs. The main 
motivation for using this approach is that the conformal techniques allow us, in principle, 
to compute the conformally extended solutions including conformal boundaries. Recall 
that in our setting, smooth conformal boundaries represent the infinite timelike future 
or past, and hence play a particular role for the cosmic no-hair picture; cf. Section 13.11 
in the first paper [5j. 

The "physical metric", by which we mean a solution of Einstein's field equations, is 
denoted by and all corresponding quantities (connection coefficients, curvature tensor 
components etc.) are marked with a tildel. The so-called "conformal metric" on the 
conformal compactification is denoted by g; all corresponding quantities are written 
without a tilde. Both metrics are related by the expression g = Q'^g, where > is a 
conformal factor. We cannot give further explanations here; some more details are listed 
in the first paper [5], and a comprehensive review is in |12j . 

We will use Friedrich's general conformal field equations in a special conformal Gauss 
gauge [m [T2I [22l El U] . In turns out that up to a rescaling of the time coordinate t of 



^Recall that, by contrast, a prime denotes a derivative with respect to the time coordinate in conformal 
time gauge in Eqs. (|2.7|l . 

''Be aware that this is the opposite notation than in [5]; there, all quantities defined with respect to 
the conformal metric are marked with a tilde. 
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this gauge, it is equivalent to a physical Gauss gauge (defined with respect to g) with 
time coordinate t. For A = 3, the rescaling has the form 



t = In 



t 



2-t 



If a smooth compact past conformal boundary J'~ exists, and if the solution extends 
to the conformal boundary in this gauge, then J'~ equals the t = 0-hypersurface where 
t — > — cxD. Under analogous conditions, J'^ is represented by the t = 2-hypersurface 
where t ^ oo. 

Now we write the evolution equations and list the unknowns. We always assume 
A = 3 in order to obtain the simplest expressions as possible. Among the unknown 
fields is a smooth frame {cj}, which is orthonormal with respect to g, and which we 
represent as follows. Due to our gauge choice, we can fix 



which is henceforth the future directed unit normal, with respect to g, of the t = const- 
hypersurfaces. Furthermore, we write 



where {ej') is a smooth 3 x 3-matrix valued function with non- vanishing determinant on 
§^ X §2. Let us define 

Wi = sin (pde + cos (j) cot 9d^, W2 = cos (f)de — sin (j) cot Od^j,, W3 = d^, (2-13) 

and from those the vector fields 



The factors 2 are chosen for later convenience. It turns out that {Va} forms a smooth 
global frame on §^ x S^. 

Having fixed the residual gauge initial data, as described in [3] , a hyperbolic reduction 
of the general conformal field equations is given by 



eo = dt 



(2.12a) 




(2.12b) 



Vi = 2(- sin cos + VFi), ¥2 = 2(sin 6 sin (/> 5p + W2) 

v^ = 2{cosedp + w^). 



(2.14) 




(2.15a) 
(2.15b) 
(2.15c) 
(2.15d) 
(2.15e) 
(2.15f) 



dtLab — —dt^Eab — XaLch, 



(2.15g) 
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for the unknowns 

b 



u 



Xab, Lab, Efe,Bfe ] . (2.15h) 



The unknowns u are the spatial components ej' of a smooth frame field {ej} as in 
Eq. (|2.12|) . the spatial frame components of the second fundamental form Xab defined 
with respect to cq, the spatial connection coefficients Vj'^, given by Tj'^Cb = Ve„ec— XacCo 
where V is the Levi-Civita covariant derivative operator of the conformal metric g, the 
spatial frame components of the Schouten tensor Lab, which is related to the Ricci tensor 
of the conformal metric by 

and the spatial frame components of the electric and magnetic parts of the rescaled 
conformal Weyl tensor Eat and Bat [121 E] , defined with respect to eo • Because the 
timelike frame field eo is hypersurface orthogonal, Xab is a symmetric tensor field. In 
order to avoid confusion, we point out that, in general, the conformal factor 17 is part of 
the unknowns in Priedrich's formulation of the CFE. However, for vacuum with arbitrary 
A and for arbitrary conformal Gauss gauges, it is possible to integrate its evolution 
equation explicitly [11], so that $7 takes the explicit form Eq. p.l5g ) for our choice of 



gauge. We note that. Eat and Bab are tracefree by definition. Hence we can get rid of 
one of the components of each tensor, for instance by substituting £'33 = —En — E22; 
we do the same for the magnetic part. The evolution equations Eqs. ()2.15ep and (|2.15fp 
of Eab and Bab are derived from the Bianchi system [12]. In our gauge, the constraint 
equations implied by the Bianchi system take the form 

D,^E^, - e'^^BdaXb'' = 0, D,^B\ + e^^E^aXb'' = 0- (2.16) 

Here, e^^^ is the totally antisymmetric symbol with e^23 — ^, and indices are shifted by 
means of the conformal metric. The other constraints of the full system above are equally 
important, but are ignored for the presentation here. Note that in Eqs. (j2.15ep . ()2.15fp 
and ()2.16p . the fields {ea} are henceforth considered as spatial differential operators, 
using Eq. ()2.12p and writing the fields as differential operators in terms of coordi- 
nates according to Eq. ()2.14p and ()2.13p . Interpreted as partial differential equations, 
these evolution equations are symmetric hyperbolic and the initial value problem is well- 
posed. Further discussions of the above evolution system and the quantities involved 
can be found in the references above. 

Friedrich's CFE allow us to use J^~^ , i.e. the t = 2-surface (or in the same way J^~), 
as the initial hypersurface. This particular initial value problem was considered in [4J. 
However, in our present application, not all solutions of interest have smooth conformal 
boundaries. In order not to exclude those solutions, we choose the t = 1-hypersurface 
as the initial hypersurface, which is a standard Cauchy surface. The hope is that this 
setup allows us to compute the complete solution including the conformal boundary if 
it exists. 

It is a standard result that there exist no globally smooth frames on x with the 
property that each frame vector field has vanishing Lie brackets with both Gowdy Killing 
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vector fields dp and 5^. The reason is given for instance in [6j. It is only possible to find 
a frame whose Lie brackets vanish for one of the two Killing vector fields, say, dp. This, 
however, has the consequence that the frame components of all tensor fields derived from 
a Gowdy invariant metric g depend on the coordinate (/> in a non-trivial manner. In order 
to reduce the evolution equations based on such a frame to 1 + 1 dimensions nevertheless, 
we can do the following [6j. It turns out to be possible to evaluate the (/)-derivative of 
every relevant unknown at, say, </> = algebraically in terms of the unknowns. Then one 
can write an evolution system which only involves ^-derivatives in space by substituting 
all ^-derivatives with these algebraic expressions. The resulting system, which we called 
1 -|- 1-system in |6j in the case of spatial S^-topology, is symmetric hyperbolic for the 
conformal field equations in our gauge. It follows from the discussion in Section El that 
a similar reduction to 1 + 1 is possible for spatial x §^-topology. The resulting system 
of equations is used exclusively in all of what follows in this paper. 

2.5 Numerical infrastructure 

It follows from Section |X] that for Gowdy symmetry, spatial §^ x and S^-topology 
have the same representation. Basically, we only need to substitute the Euler coordinates 
{Xi Pi-, P2) in the §^-case defined in Eqs. (|A.2p and (|A.4p by the coordinates {6, (p, p) in the 
§^ X S^-case, and the reference frame {Ya} defined in Eqs. (|A.3|1 by the {Va} mentioned 
before. Then, it is possible to use the same numerical technique as that which was 
worked out originally for spatial §^-topology in [6]. 

Let us repeat quickly the main ingredients of this code. By means of the Euler 
coordinates of S^, it is possible transport all geometric quantities and hence Einstein's 
field equations themselves from to T^; loosely speaking, we make "all spatial directions 
periodic". It is clear that such a map must be singular at some places. However, it is 
possible to analyze the behavior of Fourier series at the singular places and to compute 
the formally singular terms in the equations explicitly. Hence, it is not only natural 
to use Fourier based pseudospectral spatial discretization due to the periodicity in each 
spatial direction on here, but it also allows to regularize the formally singular terms 
in spectral space. This is the motivation for choosing spectral discretization in space. 
Nevertheless, a scheme to enforce "boundary conditions" [6] can be necessary in practice 
to guarantee the numerical smoothness and stability. We come back to this when we 
present our results in Section [3l 

For the time discretization, we use the method of lines. In this work, all numerical 
results were obtained with the adaptive 5th-order "embedded" Runge Kutta scheme 
from ^25j unless noted otherwise. 

2.6 Numerical computation of the initial data for the CFE 

In Section [221 we constructed initial data for the functions R, P, Q and A. Now we fix 
the initial value of the frame, i.e. the components ej' in Eq. (|2.12|) . and compute the 
corresponding initial values of u in Eq. ()2.15hp . 
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We choose the initial value of as follows. For the family of initial data constructed 
in Section 12.21 consider the frame {14} in Eq. ()2.14p and perform a Gram-Schmidt 
orthonormalization with respect to the initial conformal 3-metric. More precisely, we 
construct the matrix (e^'') from Eq. ()2.12p as an upper triangular matrix. For instance, 
this means that 63 and V3 are collinear initially. Of course there is a great freedom of 
choosing frames, and this choice is just one possibility. Note that for the following, no 
time derivative of ej^ at the initial time t = 1 needs to be prescribed. 

Now we comment on the computation of u Eq. (|2.15hp from these data. The data 
in Section 12.21 yield all spatial derivatives of the initial metric components and the first 
time derivatives. However, in order to compute u at the initial time t = 1 from these 
data, we also need second time derivatives of the data. We calculate these by imposing 
the evolution equations of Einstein's field equations at t = 1. We decided to perform all 
the computations numerically. Note that for this, the metric functions of the initial data 
in Section \2 . 21 vield formally singular terms at = 0, vr. However, we are able to compute 
these formally singular terms numerically by applying the spectral approach which was 
described above in the context of the evolution equations. In practice, we find that this 
allows us to resolve the data u with high accuracy. It turns out that machine round- 
off errors, i.e. errors introduced by the finite number representation in the computer, 
often yield the largest error contributions here. This is true in particular when the 
standard "double precision" with round-off errors of order 10^^^ on Intel processors is 
used. Hence, we decided to compute the initial data with "quad precision" of the Intel 
Fortran compiler [20j, where numbers are represented with roughly 32 digits, but which 
is software emulated and hence relatively slow. For the evolution, we switch back to 
double precision. All our numerical computations presented here have been obtained in 
this way. 

3 Results 

3.1 Choice of perturbed data 

Let us proceed by explaining our particular choices of initial data for the following 
numerical results. We present only perturbations of one generalized Nariai solution 
given by i?* = 1.0, k = 0.5 here; cf. Eq. (|2.9p . In order to perturb these Nariai data 
in an interesting manner, we choose the inhomogeneity parameter yS^^ in Eq. ()2.7p as 
non-zero, but small, and furthermore introduce a small non-zero parameter (j, by 

For the value /i = = 0, the data reduce to Nariai data according to Eq. (|2.9p . In all 
of what follows, we show the numerical results for three initial data sets given by 

/i = 0.0004667, 0.0004800, 0.0005000, 

and = 4- 10~^ in all these three cases. These choices are motivated as follows. On the 
one hand, we want to focus on "small perturbations" as a first step in order to study the 
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Fig. 1: Initial spatial dependence of H2 for the three initial data sets considered here. 

instability of the Nariai solution carefully. We believe that the choice of the parameters 
above is consistent "with this. Indeed, "we have experimented with other values of these 
parameters. In particular in a large range of values for there is no qualitative change 

in the results that follow. When we go to "very large" values T^y^ ~ 10~^, then a different 
phenomenology occurs; this interesting aspect is currently under investigation and will 
not be presented in this paper. When we, however, go to even smaller values of S^^^^ 
we get problems with numerical accuracy. Numerical errors in our runs are discussed in 
Section We remark that in this case of small perturbations, a linearization of the field 
equations around the unperturbed Nariai solution can be a reasonable approximation. 
This is currently work in progress, but will not be presented here; indeed, all numerical 
results that follow are based on the full non-linear field equations. 

In any case, smallness of the parameters is not the only motivation for our particular 
choices of data above. Recall that we expect that the instability of the Nariai solution 
can be exploited to construct cosmological black hole solutions. Our expectation is that 
this instability is controlled by the sign of the initial value of H2. For our choices of 
parameters, H2 has positive and negative parts according to Eq. (j2.1ip . see Fig. [TJ Our 
expectation for these data sets is hence that in the future, all these solutions collapse 
and form the interior of a cosmological black hole solution at those spatial places where 
H2 < initially, i.e. close to the equator of the spatial §^-factor, and expand and form 
the cosmological region with a smooth piece of at those spatial places where H2 > 
initially, i.e. at the poles of the spatial §^-factor. Changing the value of n with fixed ^ 
shifts the initial spatial profiles of H2 "vertically" in Fig. [1] leaving the shape and the 
amplitude of the curves approximately invariant. Due to this, we expect that the larger 
is, the "smaller" should be the black hole region of the resulting solution. 

3.2 The numerical results 

Now we present our numerical results based on these data sets; practical details and a 
discussion of numerical errors are given afterwards in Section 13.31 
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Fig. 2: Spatial minimum and maximum of H2 vs. time. 



Future evolution The horizontal axes in the plots in Fig. [2] represent the time co- 
ordinate t; recall that the initial hypersurface corresponds to t = 1, and t = 2 would 
correspond to the infinite timelike future. Hence, these plots show the future evolution 
of H2 for our three initial data sets. On the vertical axis, we show the maximum and 
minimum values, respectively, of H2 at any given time t. In the early phase of the evolu- 
tion, the solution behaves in accordance with our expectations. Basically, the expansion 
of the coordinate S-^-factor given by H2 becomes more and more positive where it is 
positive initially, namely at the maximum at the poles of the 2-sphere, see again Fig. [TJ 
Furthermore, it becomes more and more negative where it is negative initially, namely 
at the minimum at the equator of However, at a time t ~ 1.3, the behavior changes 
completely. The spatial profiles of H2 start to become "flatter" in the sense that the 
maximal value of H2 becomes smaller and the minimal value larger with increasing t, 
as can be seen from the figure. Eventually the solutions "make a decision" whether the 
coordinate §^-factor expands or collapses indefinitely globally in space. We give more 
evidence for this in a moment. 

We do not understand the mechanism underlying this phenomenon yet. We hope to 
be able to to shed further light on this by means of the linearization of the problem in 
future work. In any case, the numerical results suggest that there is a new instability and 
a new critical solution, in addition to the expected instability of the Nariai solutions. 
That is, there must be a critical value of fi in the interval (0.00048,0.0005). For 
/.i < /ic, the solution collapses eventually, and for /i > fic, expands globally in space. 
It would be interesting research to identify the critical solution and to study whether 
critical phenomena, which play such an important role p] for the critical collapse of 
black holes, also occur here. In any case, it is an interesting unexpected result that it 
does not seem possible to construct cosmological black hole solutions for small Gowdy 
symmetric perturbations of the Nariai solution, in contrast to the claims in [3| for the 
spherically symmetric case. 

We present further evidence for our interpretation of the numerical results now. For 
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Fig. 3: Spatial minimum and maximum the Kretschmann scalar K vs. time. 
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Fig. 4: Future late time behavior. 



this, consider the plots in Fig. [3] for the Kretschmann scalar 

where R^ippa is the Riemann tensor of the physical metric g. The curves are consistent 
with what we have just said, and confirm in particular that the collapse or expansion 
takes place globally in space eventually. The kinks in these curves can be explained as 
follows. The spatial profiles of the Kretschmann scalar in our evolutions have several 
local extrema in space which "compete" to become the global extremum. 

All the plots so far focus on the early time behavior of the solutions due to the choice 
of scales on the axes. Now let us look at Fig. HI which focuses on the evolutions at 
late times. We note that on these scales, the curves of the maxima and minima of the 
quantities are not distinguishable and hence we only show one. In the first picture, we 
show the Hubble scalar H, cf. Eqs. ()2.10p . In the eventually expanding case given by 
/X = 0.0005, we can show that the solution develops a smooth J'~^ numerically. This is 
consistent with (but not implied by) the fact that H converges to the value 1 at i = 2 
in the plot. Recall the discussion of the behavior of H at J'^ in Section [3] of the first 
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Fig. 5: Future evolution of the expansion of S^-factor. 
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paper [5] for the case A = 3. For the other two solutions, these plots confirm that they 
collapse indefinitely. This follows from a singularity theorem [2], because in both cases, 
H eventually becomes smaller than —1. The second plot in Fig. U] shows K versus t, 
and it reinforces our previous statement that the curvature of the "collapsing" solutions 
blows up everywhere in space eventually. 

Let us finish this part with a discussion of the evolution of other aspects of the 
geometry. In Fig. [5l we see the evolution of Hi, cf. Eqs. ()2.10p . i.e. the expansion of 
the coordinate §^-factor. According to this plot, we conjecture that the two collapsing 
solutions form a singularity of cigar type [29], in the same way as in the spatially homo- 
geneous case [5]. In Fig. [6l we show the maximum and minimum value of sin a, where 
the angle a is defined as follows. At a given time t = const and spatial point, a is the 
angle between the vector dp and the normal vector of the coordinate §^-factor within 
the t = const-surface. The plots suggest that a approaches zero eventually. Hence, 
loosely speaking, the Gowdy solutions become more and more polarized. Again, we do 
not understand the mechanisms underlying these curves and hope that a linearization 
will shed further light on this. Since the curves for the three solutions are almost indis- 
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tinguishable, it is a natural question whether this behavior is universal in our class of 
solutions. 



Past evolution For completeness, let us proceed with the evolution in the past time 
direction, and recall from the discussion in the first paper [S] that the unperturbed 
Nariai solutions with do < form a Cauchy horizon in the past. There is additional 
motivation from the strong cosmic censorship issue [Ij to understand what happens to 
this horizon under our perturbations. Because we consider the past time direction now, 
the initial hypersurface given by t = 1 is on the right of the following plots, and the past 
evolutions take place to the left. In the first plot of Fig. \7\ we show the maximum of 
the Hubble scalar H for our three cases of initial data together with the corresponding 
curve of the unperturbed Nariai solution. Again, corresponding curves of the minima 
are not distinguishable on these scales. We see that the four curves in the plot are 
almost the same, and hence all four cases collapse in the same way to the past. Do the 
perturbed solutions hence also develop a Cauchy horizon in the past? The second and 
third plot in Fig. [7] suggest that this is not the the curvature blows up uniformly 

for all perturbed solutions. We stress that the numerical results for the minimum of the 
Kretschmann scalar are not conclusive yet, since the value of min(i^) is still relatively 
small when the runs were stopped. Nevertheless, first signs of curvature blow up are 
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Table 1: Numerical parameters for the runs presented in Section [3.2[ 



apparent. Our observations in [1| were quite similar, and inm{K) often blew up much 
less than max(i^) close to a singularity. It is expected that this is not a geometrical 
phenomenon, but rather caused by the choice of the Gauss gauge, as we discuss there. 

3.3 Practical details about the runs and numerical errors 

Further technical details Our general numerical setup has been described in Sec- 
tion 12.51 In Table [1] now, we list more technical details about the runs in the previous 
section. The quantities Nq, Ni and /i are related to the spatial resolution. Our numerical 
runs use the simple spatial adaption technique described in [B]. After some experiments, 
the quantity X22 was chosen as the reference variable. The threshold value for the spatial 
adaption is called fi, the initial number of spatial grid points is A'o, and the number of 
spatial grid points at the stop time ti is referred to as A'^i. The following columns in 
the table describe the time discretization. We use the 5th-order "embedded" adaptive 
Runge Kutta scheme with control parameter -q. This parameter was introduced in [6] in 
order to control the desired accuracy of the time integration; the lower its value is, the 
smaller are the time steps chosen by the adaption algorithm. 

Furthermore, ho is the initial time step and hi is the time step at the stop time ti. 
In order to prevent the code from reaching unpractically small values of h, the adaption 
is switched off when h goes below hmin- One sees that for all the runs, this minimum 
value was reached eventually. Note that all numbers in the table are rounded. 

Numerical errors and convergence Prior to the numerical runs in the previous 
section, we made further tests of the code in addition to those in [6]. The choice of 
orthonormal frame in Section 12.61 has the consequence that even spatially homogeneous 
solutions "appear inhomogeneous" , in the sense that many resulting unknown tensor 
components depend on the spatial coordinates. Spatially homogeneous solutions hence 
yield a non-trivial test case for the code. These tests showed that the code is able to 
reproduce these solutions with promisingly small errors, in particular the Nariai solution 
itself. 

In order to give the reader an impression of the size of numerical errors in the 
results in the previous section, let us redo the run for the initial data set fi = 0.00048 
and ^ = 4 • 10-^ to the future with other resolutions than in Table [TJ In order to 
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study convergence more cleanly, let us switch off all adaption techniques for this. First, 
consider Fig. [8l For the same initial data, we made six runs with the spatial resolutions 
N = 100, 150, 200, 250, 300 and 350, and fixed size of time step h = 10"^. The figure 
shows the absolute values of the differences of two successive runs for the quantity e^^ 
at 6 = versus time t. We note that we have looked at other variables and seen the 
same results qualitatively. The absolute size of these differences can be interpreted as a 
measure for the size of the absolute pure numerical error for a given spatial resolution 

(or equivalently spectral truncation) on the one hand. On the other hand, since these 
errors get smaller dramatically for increasing N, we have demonstrated convergence of 
the error in our numerical results. Up to ~ 300, the numerical errors in these runs are 
hence dominated by the spectral discretization. Increasing N further, does not decrease 
the numerical error, and other types of errors become dominant, in particular the errors 
given by the time discretization and machine round-off errors. The plot also allows us 
to quantify the rate of convergence. For resolutions smaller than N ~ 250, we find that 
50 additional grid points decrease the error by a factor of approximately 3000. This 
shows that the convergence is exponential in this regime, and hence confirms that our 
numerical techniques are reliable and the numerical errors in our results in Section 13.21 
are small. The fact that spatial resolutions N ~ 300 are necessary in order to make 
spatial discretization errors smaller than other errors even at early times, when the 
solutions are very smooth in space in principle, shows that our choice of frame is not 
optimal. However, the fact that we see such a nice convergence for a quantity evaluated 
at the coordinate singularity 6 = provides particular evidence that our numerical 
regularization of the coordinate singularities mentioned in Section 12.51 works well. 

In Fig. [9l we show the same for fixed spatial resolution N = 300 and the following 
time resolutions: h = 8 ■ 10~^, 4 • 10~^, 2 • 10~^ and 1 • 10~^. As mentioned earlier, the 
errors given by the spatial discretization should be negligible for A^ = 300. Note, that 
here, instead of the adaptive 5th order Runge Kutta scheme, we use the standard 4th 
order Runge Kutta scheme. The figure confirms 4th order convergence of the errors as 
long as those are dominated by time discretization. This is the case in particular for 
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later evolution times. At very early times, however, the errors are strongly influenced 
by the machine round-off errors and so barely converge with increasing resolution. If we 
decreased h even further, the errors would be more and more dominated by round-off 
errors for longer and longer evolution times and convergence would be lost. Again, all 
this confirms that our numerical techniques are reliable and the numerical errors in our 
results in Section 13.21 are well understood and small. 

For the discussion of other important error quantities, we introduce the following 
definitions from OS]. First, we define 



with the physical Ricci tensor Rij evaluated algebraically from the conformal Schouten 
tensor Lij and derivatives of the conformal factor $7. The spatial slice at time t is 
referred to as T,t here. The indices involved in this expression are defined with respect 
to the physical orthonormal frame given by = 0,ei, and we sum over the L^-norms of 
each component. Hence, this norm yields a measure of how well the numerical solution 
satisfies Einstein's field equations Eq. (jl.ip . Second, let us define Norm^'^""''*'"-' as the 
L^-norm of the sum of the absolute values of each of the six components of the left hand 
sides of Eqs. (j2.16p at a given instant of time t. For the definition of the norm Norm^^"^, 
we again refer to [6j. The smoothness of the solution implies a certain behavior of all 
unknowns in our evolution problem at the coordinate singularities at ^ = 0, vr, and the 
quantity Norm^'^'^ is the sum of the absolute values of all quantities, that, in line with 
this behavior, should vanish at ^ = 0, vr at a given time of the evolution. 

These norms are used in the following. In addition to the pure numerical errors of the 
type discussed above, numerical relativity is plagued with the "continuum instability" 
of the constraint hypersurface in general when the constraints are propagated freely. We 
hence stress that this is not a particular problem of our numerical investigations here. 
For some evolution systems, one is able to control the constraint behavior slightly [18j, 
but a general solution to this fundamental problem has not yet been found. When the 
evolution is started with an arbitrary small violation of the constraints, then typically. 
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these violations grow exponentially, or even blow up after finite time, even for the con- 
tinuum (i.e. non-discretized) equations [13]. In Fig. [TOl we see the constraint violations 
for the runs to the future of the previous section. In addition, we show the constraint 
propagation for the case /i = 0.00048 and S^^^ = 4 • 10"^ with the same numerical pa- 
rameters as for the case /U = 0.0005 from Table [TJ In accordance with typical numerical 
runs, we see that the constraint violations grow strongly during the evolution, almost in 
the same way whether the solutions collapse or expand eventually. Increasing the ini- 
tial spatial resolution leads to higher initial constraint violations due to initially higher 
machine round-off errors. Nonetheless, the plot suggests that the constraint violations 
decrease once the discretization errors become dominant. This positive result is consis- 
tent with our observations in [3], and demonstrates that the constraint violations can 
be controlled and kept close to the continuum evolution of the constraints for arbitrary 
large evolution times as long as the errors are dominated by discretization and not by 
machine round-off errors. As discussed above, for the continuum equations, the initial 
value of the constraint violation is of fundamental importance. We have tested this for 

= 300 in the following way, as we do not show here. We have repeated some of the 
runs before with "quad precision", mentioned in Section [231 With quad precision, the 
initial size of the constraint violation is many orders of magnitude smaller than in the 
standard "double precision" case. Our numerical evolutions show that this stays true 
for the whole evolution in particular because the machine round-off errors, which would 
otherwise dominate for N = 300 and sufficiently high time resolution, are much smaller. 
Now, since the results presented in Section 13.21 are virtually unchanged when they are 
repeated with quad precision, we conclude that these results are reliable despite the ap- 
parently large violation of the constraints at late times. Concerning the violation of the 
full Einstein's field equations, hence including all constraints and evolution equations, 
the same arguments lead to similar conclusions; consider Fig. [TTl 

Let us point the attention of the reader to Fig. [T2l in order to show the order of 
magnitude of the violation of the smoothness conditions at the coordinate singularities. 
As we do not show here, these errors converge to zero, as long as the pure numerical 
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errors are not dominated by machine round-off errors. We note that none of the runs 
presented here enforce these smoothness conditions exphcitly; it is possible that this 
would improve the numerical accuracy slightly [6]. 

4 Summary and outlook 

In this paper, we studied the instability (non-genericity) of the Nariai solutions for the 
family of Gowdy perturbations. The investigations here are based on the first paper 
[5]. Our motivations to do this were two-fold. First, we were interested in the funda- 
mental question of cosmic no-hair and its dynamical realization in more general classes 
than the spatially homogeneous case considered in [5j. Second, the results of the first 
paper suggest that the understanding of the instability of the Nariai solutions in the 
spatially homogeneous case could be exploited in order to construct cosmological black 
hole solutions with in principle arbitrarily complicated combinations of black hole and 
cosmological regions. Indeed, this interesting possibility was already considered in the 
spherically symmetric case in [7], where the author claims that such constructions are 
possible. Since no non-trivial cosmological black hole solutions are known for Gowdy 
symmetry with spatial x S^-topology to our knowledge, it was our aim to address this 
open problem. 

Our results, which are obtained with the numerical technique introduced in [6j, are 
as follows. First, by making experiments with various choices of perturbations, indeed 
more than those presented in this paper, we can confirm the expected instability of 
Nariai solutions, and hence the cosmic no-hair conjecture also in the case of Gowdy 
symmetric perturbations of the Nariai solution. That is, either the solutions close to a 
Nariai spacetime collapse in a given time direction, or when they expand, they form a 
smooth conformal boundary and hence are consistent with the cosmic no-hair picture. 
Hence, our results can be seen as a generalization of the work in [5j on the one hand. 
However, of even stronger interest is that our numerical results suggest that it is not 
possible to construct cosmological black hole solutions with small Gowdy symmetric 
perturbations of the Nariai solutions. This result is unexpected, in particular it is con- 
trary to the claims for spherical symmetry. We find that the early time behavior agrees 
with the expectations. But then, the quantity H2 starts to level off and the solution 
makes a decision, whether to either expand or collapse globally in space. The underly- 
ing mechanism is not understood. Of particular interest for future research will be the 
construction and study of the critical solution and of possible critical phenomena. One 
promising approach for shedding further light on these issues is to linearize the problem, 
on the one hand around the unperturbed Nariai solution, and on the other hand around 
the hypothetical critical solution. 

Certainly, our class of initial data cannot be considered as "generic" , or to put it the 
other way around, it is not clear how "special" it is. Thus it is hard to make predictions 
for general solutions close to generalized Nariai spacetimes. We are currently working 
on a method to obtain "general" Gowdy symmetric initial data numerically. General 
Gowdy initial data would allow us to study generalized Nariai solutions in particular in 
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the standard case ctq > 0. In this hght, we understand our results here as first steps in 
an ongoing research project. 

Since we find that it does not seem to be possible to construct cosmological black hole 
solutions by means of small Gowdy symmetric perturbations of Nariai data, it is natural 
to investigate large perturbations as a next step. The hope is that the spatially local 
behavior, which is suppressed in the case of small perturbations apparently, becomes 
significant. Beyond what we have presented in this paper here, our preliminary results 
suggest that this is the case. This will be investigated in another future publication. 

We have discussed numerical errors and given some evidence that our numerical 
results are reliable. However, there is certainly room for improvements, not only in the 
numerical techniques, but also in the choice of gauge and the particular formulation of 
the field equations. For instance, in [3], we have interpreted the fact that we do not 
see spatially local behavior close to the singularities in our runs, so-called Gowdy spikes 
yy , as a reflection of the "bad" features of the Gauss gauge. That is, in this gauge, the 
solution approaches the singularity in a too inhomogeneous manner, obscuring such small 
scale structure. Hence, other gauge choices should be investigated. Another problem, 
already addressed before, is that our particular evolution system does not show optimal 
constraint propagation. Other formulations of the system should be tried and "constraint 
damping terms" should be investigated in order to improve this problem. Nevertheless, 
we have concluded above that our current numerical results can be trusted despite the 
apparently large constraint violations. 
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Appendix A Relation of the U(l) x U(l)-actions on §^ x 
and §3 

Let us consider a smooth global effective action of the group U(l) x U(l) on a 3- 
dimensional manifold. One can show, see the references in [9j, that the only compatible 
smooth compact orientable 3-manifolds are T'^, x S'^ and lens spaces. Since the 
universal cover of the lens spaces is S^, they will always be included when we speak about 
S^. In this paper, we are particularly interested in the case x However, for the 
discussion of our numerical approach, which was worked out for the §^-case originally, 
it makes sense to consider the case §^ simultaneously now. On §^ x S^, let us consider 
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coordinates {p, 6, (p) where p G (0, 2it) is the standard parameter on §^ and (9, cp) are 
standard polar coordinates on §^ 

xi = sin ^ cos = sin sin xs = cos9. i-^-^) 

In writing these coordinate expressions, we assume that is embedded in the standard 
way into with Cartesian coordinates (xi, X2, X3). On S^, we consider Euler coordinates 
(x, Ai, A2) with the same conventions as in [4J, namely 

xi = cos — cos Ai , X2 = COS — sin Ai , 

X X (^-2) 

X3 = sin — cos A2, X4 = sin — sin A2, 

where we assume the standard embedding of into similar to the above. Here, 
X £ (0, vr) and Ai,A2 G (0,27r). In terms of these coordinates on §^ x and 
respectively, we can write a representation of the action of the group G = U(l) x U(l). 
For X S^, one has 

^' : G X (§1 X §2) ^ §1 X ((ni, U2), (p, 6, 0)) ^ (p + ui, 6*, + U2) 

for (^1,^2) G U(l) X U(l) and (p,9,(j)) £ x S^, writing this, we always assume the 
standard identification of the groups U(l) and Hence a basis of generators of the 
action are the coordinate fields dp and d^. These are globally smooth vector fields on 
X §2 and one can check that this action is global, smooth and effective. The action 
degenerates at those points where the vector field has a zero, namely at the poles of 
the S^.factor given by = 0, vr. In the case of S^, we have the following action 

^ : G X §3 ^ {{ui,U2), ix, Ai, A2)) ^ (x, Ai + ui, A2 + U2). 

The generators of the group are the coordinate fields dx^ and dx^ . Here, as well, these are 
globally smooth vector fields on S'^. Indeed, the action is global, smooth and effective. 
The action degenerates where either dx^ or 8x2 have a zero, which is the case at x = 0, vr. 
Now, it is a fact, quoted in [9] , that all other smooth effective global actions of U(l) x U(l) 
on any of these manifolds must be equivalent, i.e. can only differ by an automorphism 
of the group or by a diffeomorphism of the manifold to itself. Hence it is sufficient to 
have a single representation of the action. 

We will now formulate the actions above in an equivalent, but more geometrical 
manner. We note that the fields Wi, W2, W3 defined in Eqs. (j2.13p . together with 

6 := dp, 

form a basis of the Killing algebra in the spatially homogeneous case on xS^. Moreover, 
the action of the group U(l) x U(l) is generated by {^3,^3}, i.e. the translation vector 
field along the S^-factor of the manifold and a rotation of the S^.factor of the manifold 
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corresponding to an element of the Lie algebra of SO (3). In the S'^-case, let us introduce 
the vector fields 





= 2 sin + 2 cos pi (cot x dp^ — 


esc X Spa) , 


(A.3a) 


Y2 


= 2 cos pid^ — 2 sin pi (cot x dpi — 


cscx^pj , 


(A.Sb) 




= 25,,, 




(A.3c) 


Zi 


= —2 sin p2d-^ — 2 cos p2 (cot x Qp\ 


- CSC X Spa) , 


(A.Sd) 


Z2 


= 2 cos p2d^ — 2 sin p2 (cot x dp^ — 


CSC X Spa) ' 


(A.3e) 


Zz 


= 25p,, 




(A.3f) 



in terms of the Euler angle coordinates above, but with pi and p2 defined by 

Ai = (pi + P2)/2, A2 = {pi - P2)/2. (A.4) 

The definitions of these fields are in agreement with the conventions in [6l 0] . These vec- 
tor fields have the property that they are globally smooth on S'^ and are invariant under 
the standard left- and right-actions, respectively, of the group SU(2) on S^. Furthermore, 
they satisfy 

3 3 
[Ya, n] = 2 ^ iiabcYc, \Za. ^fe] = 2 r/,fe,Z„ Zb] = 0. 

c=l c=l 

Here r]ahc is the totally antisymmetric symbol with 77123 = 1. Note that the collections 
{Yi,l2 5^3} and {^1,^2,^3} are global frames on S^. A basis of the generators of 
the U(l) X U(l)-action on S'^ is {13, Z3}. We remark that both ^3 and Z3 are nowhere 
vanishing vector fields. They become, however, collinear at x = 0, vr, and the U(l) xU(l)- 
action degenerates. 

Now, recall that both x and are principal fiber bundles over with structure 
group U(l). In the x §^-case, it is the trivial bundle with bundle map 

$i:S^xS2^S^ ip,q)^q. 

In the S^-case, it is the Hopf bundle and the bundle map can be given the representation 

^>2 : ^ S^, (xi,x2,x3,x4) ^ (yi,y2,y3) 

= (-2(-XiX3 -I- X2X4), 2(X2X3 -I- X1X4), xf -I- X2 - x| - X4). 

Here the notation is such that (xi, X2, X3, X4) G S'^ C with standard Cartesian co- 
ordinates on M^, and (yi, 2/2) 2/3) £ C M^. When we write this map in terms of the 
Euler angle coordinates for S^, given by Eqs. ()A.2p and ()A.4p . and polar coordinates on 
given by Eq. (|A.ip (with Xj substituted by y^), then the Hopf map has the simple 
representation 

^2 : {x,Pi,P2) ^ {0,(t)) = {x,Pi)- 
In order to simplify the notation we write now 
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1. N for §1 X §2 or 

2. r] for ^3 or Z3, and a for 2W3 or Y3, 

3. ^ for or ^2, 

in the following. Furthermore, we write G = U(l) x U(l), and Gi = U(l) to denote 
the subgroup of G generated by r/. Recall that r] is nowhere vanishing and note that its 
integral curves are closed circles. We notice that in both cases, the field r] is tangent to 
the fibers of the bundles, and hence Gi becomes the structure group, so that N/Gi = S^. 
Note that, in particular, $*(T is a smooth global vector field on S^, and indeed equals the 
coordinate vector field 25,^. This fact holds irrespective of the choice of either x 
or S^, and thus clearly demonstrates how similar the two cases are at this level. 

The bundle maps $ can be used as follows. Assume that we want to solve a set 
of partial differential equations on N, possibly with an additional time function t. We 
suppose that all unknowns and coefficients are smooth functions on constant along ij, 
and that all differential operators in the equations originate in smooth vector fields on N 
whose Lie brackets with r] vanish. Then it turns out that there is an equivalent system of 
equations on as follows. Each such unknown and coefficient can be identified uniquely 
with a smooth function on by means of the bundle map <I>. Furthermore, each such 
vector field yields a unique smooth vector field on (which certainly has zeroes) when 
pushed forward with Most importantly, the solution of these new equations on S^, 
if it exists, hence yields a unique corresponding solution of the original problem on N. 
Thus without loss of information, we are allowed to "transport such geometric problems 
from N to S^" along the bundle map 
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